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' Abstract 

The cross section for the inclusive photoproduction of large-pi 1 D*^ mesons is 
calculated at next-to-leading order, adopting different approaches to describe the 
fragmentation of charm quarks into D*^ mesons. We treat the charm quark accord- 
ing to the massless factorization scheme, where it is assumed to be one of the active 
flavours inside the proton and the photon. We present inclusive single-particle dis- 
tributions in transverse momentum and rapidity, including the contributions due to 
both direct and resolved photons. We compare and assess the various implementa- 
tions of fragmentation. We argue that, in the high-p^ regime, a particularly realistic 
description can be obtained by convoluting the Altarelli-Parisi-evolved fragmenta- 
tion functions of Peterson et al. with the hard-scattering cross sections of massless 
partons where the factorization of the collinear singularities associated with final- 
state charm quarks is converted to the massive-charm scheme. The predictions thus 
obtained agree well with recent experimental data by the HI and ZEUS Collabora- 
tions at DESY HERA. 
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1 Introduction 



Heavy-quark production in the ep colliding-beam experiments with DESY HERA offers 
a novel way of testing perturbative quantum chromodynamics (QCD). First results for 
the charm-quark photoproduction cross section a(jp — > cc + X) were presented by ZEUS 
and HI 0, and compared with a next-to-leading-order (NLO) calculation ||. In this 
calculation, the massive-charm scheme has been adopted, in which the charm-quark mass 
m ^> Aqcd acts as a cutoff and sets the scale for the perturbative calculation. The cross 
section factorizes into a partonic hard-scattering cross section multiplied by light-quark 
and gluon densities [[|. In this factorization approach, the only quarks inside the proton 
and the photon are the light ones. Thus, in the massive-charm scheme, the number of 
active flavours in the initial state is rif = 3, while the massive charm quark appears only 
in the final state. For the prediction of the total charm photoproduction cross section, for 
which experimental results from HERA were presented in ]I]j2| , this is the only possibility. 
Actually, m ~ 1.5 GeV is not very large compared to Aqcd, so that the validity of this 
approach is not obvious. 

With the advent of new measurements, by HI |2| and ZEUS ||, of the differential 
cross section d 2 a / ' dy dp? of inclusive D*^ production, where y and px are the rapidity 
and transverse momentum of the D*^ mesons, respectively, we have the possibility to 
test the theory in a different regime of scales. The experimental differential cross sections 
extend up to pr = 12 GeV, so that, in contrast to the total-cross-section calculations, pt 
rather than m should be considered as the large scale. Then, in NLO, terms proportional 
to a s ln(p T /m 2 ) arise from collinear gluon emission by charm quarks or from almost 
collinear branching of gluons or photons into charm-anticharm pairs. For large enough 
Pt, these terms are bound to spoil the convergence of the perturbative series and cause 
large scale dependences of the NLO result at pt 3> m. The proper procedure in the regime 
Pt 3> m is to absorb the terms proportional to a s \n(p T /m 2 ) into the charm distribution 
functions of the incoming photon and proton and into the fragmentation functions (FF's) 
of charm quarks into charmed hadrons. Of course, to perform this absorption, one needs 
information on the charm contribution in the parton density functions (PDF's) and FF's. 

An alternative way of making reliable predictions at large pt is to treat the charm 
quarks as massless partons. The collinear singularities corresponding to the a s ln(p T /m 2 ) 
terms of the massive-charm scheme are then absorbed into the charm-quark PDF's and 
FF's in the same way as for the lighter u, d and s quarks. This massless approach was 
proposed in || and first applied to the production of large-pr hadrons containing bottom 
quarks in pp collisions by Lampe and by Cacciari and Greco ||. Subsequently, it was 
employed in two independent studies of charm-quark photoproduction HJT0| . In these two 
investigations, the y and pt distributions in the massless and massive approaches were 
compared with each other, making different assumptions concerning the initial state as 
well as the FF's in the massless approach. In ||, low-Q 2 electroproduction was considered, 
while [10| was concerned with photoproduction with fixed photon energy. In |J, the FF 



of the charm quark into charmed hadrons was approximated by 5(1 — z), where z = Pd/Pc 
is the scaled momentum of the charmed hadron D, whereas the authors of pi employed 
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the perturbative FF's (PFF's) of fTT |, which they evolved from the starting scale m to 



the appropriate higher scales of order according to the usual Altarelli-Parisi evolution 
equations. 

It is the purpose of this work to remove the restriction to a scale-independent 5-type 
FF of the charm quark made in our previous study with Kramer |9[] and to adopt more 
realistic descriptions of charm-quark fragmentation, including evolution to higher scales. 
Having extensively compared the massless and massive approaches in [[|, we shall now 
focus our attention on the massless approach, which is likely to be more reliable in the 
large-pr range, in which we are primarily interested here. Specifically, we shall consider the 
following models of charm-quark fragmentation: (i) 5-function-type FF without evolution 
|| for reference; (ii) PFF Jll| with and without evolution; and (iii) Peterson fragmentation 



I2|| with evolution. Choice (iii) will be considered as the most realistic one and will be 



used for comparisons with the recent HI and ZEUS || data. 

The outline of our work is as follows. In Section 2, we shall shortly describe the 
basic formalism of charm-quark fragmentation, discuss the transition from massless to 
massive factorization of final-state collinear singularities, and specify our assumptions 
concerning choices (i), (ii) and (iii). In Section 3, we shall investigate numerically the 
differences between the various implementations of charm-quark fragmentation discussed 
in Section 2. There, we also confront our final predictions based on the Peterson FF with 
the HI and ZEUS data. Our conclusions will be summarized in Section 4. 



2 Fragmentation function approach 

In this section, we describe the underlying assumptions for the massless approach, making 
various choices for the FF's of the (anti)charm quark into D*^ mesons. As is well known, 
two mechanisms contribute to the photoproduction of charm quarks in ep collisions: (i) In 
the direct photoproduction mechanism, the photon couples directly to the quarks, which, 
besides the massless u, d and s quarks, also include the massless c quark. In this case, 
no spectator particles travel along the momentum direction of the photon, (ii) In the 
resolved photoproduction mechanism, the photon splits up into a flux of u, d, s, c quarks 
and gluons, which then interact with the partons coming from the proton leading to the 
production of charm quarks at large p?- The contributing parton- level processes are the 
same as in the case of charm-quark production in hadron-hadron collisions. The charm 
quark is accompanied by a spectator jet in the photon direction. Therefore, the 7^ cross 
section depends not only on the PDF's of the proton but also on those of the photon. The 
main difference relative to the usually considered massive-charm scheme is that the charm 
quark also contributes via the PDF's of the proton and photon, i.e. charm is already an 
active flavour in the initial state. In other words, there is no essential difference between 
the treatment of the truly light quarks u, d, s and the charm quark in the initial state. 
This approach is justified if a large scale is governing the production process. In our case, 
this is the Pt of the D*^ mesons, with px 3> m. In this region, non-singular mass terms 
are suppressed in the cross sections by powers of m/pT- The important mass terms appear 
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if the virtuality of the charm quark is small. This happens in the initial state if the charm 
quark is emitted from the proton or photon, and in the final state if the partons emit a 
cc pair. These two contributions lead to the a s ln(p^/m 2 ) terms in the massive scheme. 
In our massless approach, they are summed into the scale-dependent PDF's and FF's of 
charm quarks into D*^ mesons. 

Thus, when we proceed to NLO, the following steps are taken: (i) The hard-scattering 
cross section for the direct- and resolved-photon processes are calculated in the massless 
approximation with rif = 4 active flavours. The collinear singularities are subtracted ac- 
cording to the MS scheme. Since the charm quark is taken to be massless, the singularities 
from its splittings are subtracted as well. 

(ii) The charm quark is accommodated in the PDF's of the proton and photon as a 
light flavour. The finite mass of the charm quark is taken into account by including it in 
the evolution in such a way that its PDF's are only non-vanishing above a scale set by its 
mass. 

(iii) The FF's characterize the hadronization of the massless partons, including the 
charm quark, into mesonic or baryonic states containing the charm quark. In our case, 
these are the D*^ mesons. Similarly to the fragmentation into light mesons, these FF's are 
basically non-perturbative input and must be determined by experiment, for example from 
the cross section of inclusive D*^ production in e + e~ annihilation [13|. This information 
determines the FF's at some starting scale, /i . 

An alternative that exploits the fact that the charm mass satisfies m 3> Aqcd is the 
calculation of universal starting conditions for the FF's within perturbative QCD at a 
scale no of order m. The FF's thus obtained are the PFF's mentioned above. In |TI 
these starting conditions were calculated at NLO in the MS scheme. They read 



D? = D?(x,fio) = 5(1 -x) + ^4r^C F f ' 



In - 2 ln(l - x) - 1 

m 2 



2tt [1-x 
—T f [x +(l-x)]ln- 
D?^x,to)=I% A Jx,to) = 0. (1) 



Df(x, ho) = Df(x, /i ) = ^^T f [x 2 + (1 - xf] In ■ , . 



Here, Df (Df) refers to the fragmentation of parton i into a hadron D (D) containing a 
c (c) quark, Cp = 4/3, Tf = 1/2, and q — u,d, s. The normalization of these FF's is such 
that D stands for the sum of all hadrons with a c quark. To obtain the FF into a specific 
charmed meson, for example the D* + , we must include the branching ratio B(c —>■ D* + ). 

(iv) For the higher scales fi > Ho-, the PDF's, the non-perturbative FF's and the PFF's 
are evolved in NLO up to the chosen factorization scale (which we take to be of order 
Pt) via the Altarelli-Parisi equations and convoluted with the NLO hard-scattering cross 
sections. The charm-quark mass m then only enters in terms of the starting conditions 
of the PDF's and FF's. 

As already mentioned, the large logarithmic terms proportional to ln(p^/m 2 ), which 
appear in the massive scheme, are resummed in this approach. To be specific, they are 
effectively split into two parts. One part proportional to ln(p^//x 2 ), where \x is some generic 
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factorization scale, appears in the hard-scattering cross sections having no dependence on 
m. This part may be kept small by choosing p pa p?. The other part, proportional 
to ln(/i 2 /m 2 ), is absorbed into the PDF's and FF's. The logarithm ln(p 2 /pl), which is 
large if po ~ m and p ~ pr, is incorporated via the evolution equations and therefore 
resummed. The residual term proportional to ln(pl/m 2 ), connected with the starting 
conditions in ([[]), is treated with fixed-order perturbation theory in the case of the PFF's, 
or is part of the non-perturbative input in the case of general FF's. 

The PFF's given in terms of the starting conditions ([!]) suggest some simplifications 
which we consider as reasonable approximations for the PFF's. The dominant term in (|I|) 
is the 5 function expressing the fact that, due torn > Aqcd, the charm quark transforms 
into the charmed hadron having the same momentum, any binding effect of the light 
quark of order Aqcd being neglected. The approximation D®(x,p) = 5(1 — x), denoted 
A for later use, was considered in our earlier work ||. This approximation neglects the 
important ln(p 2 /p 2 ) ) terms induced in the FF's through the evolution. Unless otherwise 

stated, we shall choose po = 2m and p = 2mr, where mx = \j m 2 + p T is the transverse 
mass of the produced hadron. 

Approximation B consists of keeping all terms in (H), but with the replacement po = p. 
This approximation includes the large ln(/x 2 /m 2 ) terms to fixed order in a s , i.e. not in 
the resummed form because the FF's are not evolved up to the scale p via the evolution 
equations. Except for the fact that the proton and photon PDF's still have four active 
flavours, this approximation is, from the conceptual point of view, most closely related 
to the massive scheme 0,0], in the sense that the important mass terms of the massive 
calculation, namely the logarithmic ones related to the collinear would-be singularities, 
are properly included, while the non-singular power terms in m/pT are omitted. We stress 
that the latter are indispensible if px is of order m, and must be included if total cross 
sections are to be calculated; approximation B is then bound to break down. However, 
our goal is to improve the theoretical description at px 3> m. Since, in case B, the 
FF's are not evolved to higher scales according to the evolution equations, the terms 
proportional to a s in (|I|) may be shifted to the hard-scattering cross section, which we 
actually choose to do. In fact, case B may be viewed as a change of factorization scheme, 
in the sense that the factorization of the final-state collinear singularities associated with 
the charm quark is adjusted so as to match the finite- m calculation. This is an equivalent 



interpretation of the matching procedure proposed in ||TTJ between the massless- charm 
calculation in connection with the PFF's and the massive-charm calculation without FF's. 
Specifically, we substitute in the hard-scattering cross sections a s (p)P < ^' T \x) \n(s/p 2 ) — > 
a s (p)Pci' T \x) ln(s//i 2 ) —a s (p )d ci (x), where P^' T \x) are the timelike LO % — >• c splitting 
functions |L4| and d ci (x) may be gleaned from ([!]): 



d r Jx) = ~Cf 



l+x 2 



x 



111 



m 2 



-21n(l-x) - 1 



p(*/)( x )ln^- + C F {-26(l-x) + 2 



x 



+ 4 



ln(l — x) 



x 
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- (l + z)[l + 21n(l-x)] 

d cg (x) = -T f [x 2 + (1 - x) 2 ] In 4 

m 2 

Mo 



dcg(a;) = 4g(^) = 4c(a;) = 0, (2) 



with q — u,d, s. The d 5 i functions emerge by charge conjugation. Since the PFF's of ([[]) 
refer to the MS factorization scheme for the collinear singularities associated with outgoing 
massless- charm lines, so do the dij functions of (Q). In case B, hi(p^/m 2 ) ~ ln(p^/m 2 ) 
is large. Leaving aside the terms in (|j) that are not enhanced by this logarithm, case 
B obviously amounts to choosing the final-state factorization scale in the massless- charm 
hard-scattering cross sections to be /i = m and omitting the d^ functions, i.e. to choosing 
the final-state factorization scale in case A to be \x = m. In our implementation of case 
B, we have D^(x,p) = 5(1 — x) for all factorization scales p, just as in case A. 

Obviously, case B may be refined by resumming the large logarithms lnQo^/m 2 ) via 
evolution. To this end, in case C, we take D®(x, p ) = 5(1 — x) as the input distribution 
at the starting scale po for the NLO evolution up to p. At the same time, we keep the dis- 
functions, with the large logarithm stripped off, i.e. with p$ now of order m rather than 
Pt, in the hard-scattering cross sections. 

An alternative refinement of case B may be obtained by taking the full PFF distribu- 
tions of ([!]) as input for the evolution up to scale p, and convoluting the outcome with 
the massless-charm hard-scattering cross sections (without d^ functions). This approach, 
which we denote by D, corresponds to the PFF approach advocated in |]TTJ for inclusive 
heavy-quark production via e + e~ annihilation. It was applied to inclusive bottom-quark 
production in pp collisions in || and to inclusive charm-quark production in 7p and 77 
collisions in |T0|JT5[| , respectively. 



However, it is not at all clear that the PFF's give the correct description for the 
fragmentation of the charm quark into charmed hadrons, since the charm quark is only 
moderately heavy. Although the average value (x) near p = /i is not precisely known, 
it is bound to be (x) < 1. Most measurements of (x) are at larger scales. In the case 



of c -> D*+, we have (x) = 0.495 +J;J}J(stat.) ± 0.007(sys.) from ALEPH |TJ and (x) 



0.515+J;SJJ(stat.) ± 0.010(sys.) from OPAL 0at/i = 91.2 GeV, and (a?) = 0.73 ± 0.05 
at p, — 10 GeV [[18], indicating that (x) < 1 at /x = po as well. Therefore, the distribution 
by Peterson et al. is usually considered to be a better approximation for the FF at 



the starting scale po- It has the form 



x(l — x) 



£>?(*, /io) = AT r " V \ 2 7 i2 - ( 3 ) 
[(1 — x) 2 + ex\ 

Apart from the normalization A/", it depends just on the parameter e, which is related to 
(x). Note that (|3p turns into a 5 function in the limit e — > 0. We choose e = 0.06, which 
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is the central value obtained in [[19] for charm quarks. In case E, we adopt (|3|) with 
A r-i_ (e 2 ~6e + 4) 




arctan = + arctan = > H — In e + , (4) 



(4-e)V4e-e 2 [ V4e - e 2 

so that Jq dx Df?(x, fio) = 1. For e = 0.06, one has Af ~ 0.625. This simplifies the 
comparison with the fragmentation models A-D. In Section 3.4, where we compare our 
predictions with the HERA data, we must adjust the normalization of our calculation by 
including the measured branching fraction B(c — > D* + ) fl7| . 

The Peterson FF for c — > D* + fragmentation of (§) has also been used in |20] for 



calculating the production of D*^ mesons in 7p collisions based on the purely massive 
scheme. In our approach, we need the evolved FF's at the scale /x, with (||) describing 
the input at the starting scale /x . We solve the evolution equations using moments in 
Mellin space. For the reader's convenience, the moments of (|3]) are listed in Appendix A. 
They are expressed in terms of hypergeometric functions of complex arguments, which 
are conveniently calculated from appropriate series expansions, also written down in Ap- 
pendix A. To account for the fact that the Peterson FF describes the transition from a 
massive charm quark into a physical D* + meson, whereas, for reasons explained in the In- 
troduction, we work with hard-scattering cross sections where the charm quark is treated 
as a massless flavour, we incorporate the d^ functions of (0) in our analysis just as we do 
in case C. This means that case E emerges from case C by replacing the 5-function-type 
input distribution at scale /x with the Peterson formula of (^|). 

An alternative approach would be to include the terms relating the massless and 
massive schemes in the evolution of the FF as we do in version D, where we evolve the full 
PFF's of (p. In other words, one would convolute the PFF's with the Peterson formula (|3]) 
at the scale /io, evolve these modified initial distributions at NLO up to the scale /x and 
employ the FF's thus obtained in the usual massless- charm MS analysis (without dij 
functions). In the following, this possible alternative for scheme E will be called scheme 
F. Both of these schemes properly resum the large logarithmic terms of the type ln(/x 2 / /x 2 ,) 
which emerge from (H) as /xq is driven up to /x. In scheme E, just the latter logarithms get 
resummed, while the non-logarithmic terms of (|1|) are treated in fixed-order perturbation 
theory. By contrast, in scheme F, also these non- logarithmic terms are included in the 
evolution. Scheme E is minimal in the sense that is allows us to include in the evolution 
just the long-distance part of the c — ► D* + fragmentation process, which is encoded in 
the Peterson ansatz (|3|). We shall adopt scheme E whenever we confront theoretical 
predictions with HERA data. In the limit e — > 0, cases E and F coincide with cases C 
and D, respectively. For a given value of e, one expects the difference between schemes E 
and F to be similar to the one between schemes C and D. Since schemes C and D do not 
yield physical observables, their difference is phenomenologically inconsequential. By the 
same token, it is unreasonable to expect that a cross section evaluated in schemes C or D 
should coincide with the corresponding charm-quark production cross section obtained in 
the all- massive scheme. In the case of D* ± photoproduction, this is even more apparent, 
since the cross sections of schemes C and D exhibit a significant sensitivity to the charm 
content of the resolved photon (intrinsic charm), which does not enter the all- massive 
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calculation at all. On the other hand, the difference between schemes E and F is to be 
compensated by a corresponding shift in the e parameter of @, i.e. e is scheme dependent 
at NLO. For instance, if one wishes to use the value of e determined through a NLO fit to 
e + e _ data as input for a NLO prediction of photoproduction, one has to make sure that 
both analyses are carried out in the same scheme. Only then, one may obtain a coherent 
description of high-energy e + e~ data and high-p^ ep data. We emphasize that, on the 
basis of an all-massive calculation including fragmentation, this cannot be achieved in any 
meaningful way, due to the presence of sizeable logarithms at fixed order. A consistent 
procedure within the all-massive framework would be to determine e from a NLO fit to 
low-energy e + e~ data on inclusive D*^ production, which has not been done so far. 

A second alternative to case E has been adopted in |21] to describe charmed-meson pro- 
duction at the Tevatron. In plfl , the evolved PFF's are convoluted with some unevolved x 
distribution of the form N(l — x) a x^ . By contrast, we assume that the non-perturbative 
x distribution of the D* + mesons, where x is the longitudinal-momentum fraction w.r.t. 
the quasi-massive parent charm quarks, is subject to evolution. 

The calculation of the hard-scattering cross sections proceeds along the lines of our 
previous works |||22| on the basis of the direct- and resolved-photon hard-scattering cross 



section obtained in [23] and [24|, respectively. 



3 Results 

This section consists of four parts. First, we specify our assumptions concerning the 
proton and photon PDF's as well as the equivalent photon approximation. Second, we 
describe the various degrees of refinement in the treatment of fragmentation. Third, we 
discuss various theoretical uncertainties related to the proton and photon PDF's and the 
dependence on scales and parameters intrinsic to the FF's. Finally, in the fourth part, we 
present our predictions for the inclusive D* ± photoproduction under realistic kinematic 
conditions corresponding to the HI and ZEUS experiments. 



3.1 Input information 

For the calculation of the cross section d 2 a / 'dyi^dp^, we adopt the present standard 
HERA conditions, where E p = 820 GeV protons collide with E e = 27.5 GeV positrons in 
the laboratory frame. We sum over D* + and D*~ mesons. We take yi a b, which refers to 
the laboratory frame, to be positive in the proton flight direction. The quasi-real-photon 
spectrum is described in the Weizsacker- Williams approximation with the formula 



/ 7 / e (x) ^ 



1 + (1 - x) 2 Q 2 max 2 / 1 
i '- In + 2m 2 e x — 




(5) 



where x = E y /E e , Q 2 min = m 2 x 2 /{l — x), a is Sommerfeld's fine-structure constant and m e 
is the electron mass. In the case of ZEUS, where the final-state positron is not detected 
(ZEUS untagged), one has Q 2 max = 4 GeV 2 and 0.147 < x < 0.869, which corresponds 
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to •yp centre-of-mass (CM) energies of 115 GeV < W < 280 GeV ||. We shall use these 
conditions for the evaluation of the cross sections in the next two subsections. Besides the 
untagged case, the HI Collaboration has also presented data where the positron is tagged. 
The corresponding parameters will be specified later. Our default sets for the proton and 



photon PDF's are CTEQ4M 0, with = 296 MeV, and GRV |2§. We work at NLO 
in the MS scheme with rif = 4 flavours. We identify the factorization scales associated with 
the proton, photon and final-state hadron, M p , M 7 and M h , respectively, and collectively 
denote them by Mf. Our standard choice of renormalization and factorization scales is 
/i = ^rriT and Mf = 2£m^ with £ = 1. If the FF's are evolved, we take the starting scale 
to be fiQ = 2m with m — 1.5 GeV. We calculate a s (fi) from the two- loop formula with 
equal to the value used in the proton PDF's. 



3.2 Refinements in fragmentation 

In our previous work U, we made the most simple ansatz for the fragmentation of charm 
quarks into charmed hadrons, namely D®(x, fi) = 5(1— x) for all scales /i. In Section 2, this 
was referred to as case A. In Figs. |l] and || we show the predictions for da/dpT integrated 
over —1.5 < yi a b < 1 and d 2 a j dy^dpx at pt = 7 GeV, respectively, according to this 
case together with the refinements B-F discussed in the same section, adopting the ZEUS 
untagged conditions. Cases B and C differ from case A in that they take the evolution into 
account. In the rough approximation of case B, the evolution is incorporated by choosing 
/io = Mh in the definitions (0) of the d^ functions, which correct for massive final-state 
factorization in the hard-scattering cross sections. This generates large logarithms of the 
form ln(M|/m 2 ), which are resummed in case C by properly evolving the 5-type FF's 
from the starting scale /iq = 2m up to Mh- In case C, the transition to the massive 
factorization is performed by adding the dij functions with // = 2m. In contrast to 
the case B, case C includes also contributions from gluon fragmentation and, to a lesser 
extent, from light-quark fragmentation. In case D, the hard-scattering cross sections in 
the massless factorization scheme (d^ = 0) are convoluted with the evolved PFF's of ([!]). 
The effect of the d^ functions is then contained in the PFF's. 

We see in Fig. [I] that, at low px, the latter three cases (B, C and D) give very similar 
cross sections da/dpx, whereas there are appreciable differences at large p?, especially 
between cases C and D. Case D exhibits the strongest fall-off. Here not only the S 
function, but also the finite, perturbative terms of (|l]) are included in the evolution. Case 
B, where evolution is accounted for only in the leading-logarithmic approximation, has 
the smallest decrease of da/dpx with increasing px- Case C, where only the 5 function 
is evolved, lies between cases B and D, as one would expect, but closer to case B. We 
emphasize that the evolution of the FF's is performed up to NLO (cases C and D). We 
conclude from this study that the evolution of the FF must be fully taken into account 
in order to obtain realistic predictions for the px distribution. The unevolved 5-function 
fragmentation considered in (case A) is unreliable for large px- 

In 0, case A was compared with the massive NLO calculation, adopting the previous 
ZEUS conditions, which are very similar to the present ones used here. The massive NLO 
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result turned out to be roughly 45% of the case- A result. On the other hand, from Fig. |1] 
we see that, at px = 3 GeV (12 GeV), the case-f? result amounts to 67% (54%) of the 
case- A result. In other words, at pt = 3 GeV (12 GeV), the massive NLO result is roughly 
2/3 (5/6) of the case-I? result. This may be partly attributed to the negative contribution 
from the power terms in m/px, which are not included in case B. The magnitude of these 
terms rapidly decreases with increasing px- On the other hand, in contrast to the massive 
calculation, case B receives contributions from the charm quarks intrinsic to the proton 
and photon. These contributions will be investigated in Figs. |] and |5], respectively. We 
shall see that the charm content of the photon is particularly important. 

The rapidity distributions at pt = 7 GeV for the cases A, B, C and D are shown in 
Fig. 0. Here the pattern is somewhat different. In the y^ region where the cross section 
d 2 a / dyi ab dpT is maximal, i.e. —1.5 < yi a t < 1, the distributions of cases A, B, C and D 
have very similar shapes. The overall normalizations just differ according to the hierarchy 
in Fig. 0. For the larger rapidities yi a b > 2, the cross sections for cases A, B, C and D 
almost coincide, i.e. here the evolution of the FF's has little effect. 

In Figs. [I] and 0, we have also included the predictions based on the Peterson FF 
evolved to larger scales in NLO, corresponding to cases E and F. This FF is also nor- 
malized to unity, as all the other FF's, in order to facilitate the comparison. That is, 
in the limit e — > 0, the results of cases E and F coincide with those of cases C and D, 
respectively. As anticipated in Section 2, for finite e, the relative difference between cases 
E and F is very similar to that between cases C and D. In particular, the px distributions 
for cases E and F in Fig. [I] converge towards the \ow-px end, where the non-logarithmic 
terms of the PFF's as implemented in case F are just mildly affected by the evolution. 
The main effect of the Peterson FF with finite e is to reduce the overall normalization of 
the respective calculations. 

In the following investigations, we shall stick to case E and employ the NLO-evolved 
Peterson FF. This description of charm fragmentation is more realistic than the PFF's 
because the charm quark is only moderately heavy. Compared with the massive-charm 
scheme, the massless-charm scheme shows a completely different decomposition of the 
cross section into direct- and resolved-photon contributions. This was analysed in detail in 
our earlier work ||. Whereas in the massive approach direct photoproduction is dominant 
at large px, in our massless scheme both contributions are of the same order of magnitude 
0. This is demonstrated in Fig. |3] for d 2 a / dyi^dpx at px = 7 GeV. Both cross sections 
peak approximately at = 0. The resolved cross section has in addition a shoulder at 
larger y^, which originates from the gluon part of the photon PDF's. The peak at yi a b = 
is due to the charm part, as may be inferred from Fig. [|. Of course, we must bear in 
mind that the decomposition of the photoproduction cross section in direct- and resolved- 
photon contributions is ambiguous at NLO; it depends on the factorization scheme and 
scale, which we take to be the MS scheme and M 7 = 2mx, respectively. However, the 
sum of the two contributions is unambiguously defined. 
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3.3 Theoretical uncertainties 



Before we compare our results with recent experimental data from HI and ZEUS, we 
investigate several theoretical uncertainties, which might be relevant. First, we consider 
the influence of the proton PDF's. In Fig. [|, the cross section d 2 a/dyi a b dpx at pr = 7 GeV 
is plotted for our standard CTEQ4M proton PDF's and compared with the result 



obtained using the recent proton PDF set MRS(G) by the Durham group 127], with 



A^L = 254 MeV. The cross sections for these two choices are almost identical. From this 
we conclude that our predictions are insignificantly influenced through our choice of the 
proton PDF's. In Fig. [|, we also show the influence of the charm distribution inside the 
proton. If it is taken away from the CTEQ4M PDF's, the cross section hardly changes. 
Only for larger y^, where the cross section is small, do we see the effect of the charm 
content of the proton. 

A similar study w.r.t. the photon PDF's is presented in Fig. ||. Here we compare 
d 2 a / dyi ab dp T at = 7 GeV, evaluated with our standard GRV photon PDF's, with the 
calculation based on the recently published new version of the photon PDF's by Gordon 



and Storrow f28j. The cross sections differ somewhat, by typically 10-20%, over the whole 



yiab range. This means that very accurate data are needed if the differences between the 
GRV and GS photon PDF's are to be disentangled. In Fig. [|, we have also plotted the 
GRV prediction with the charm distribution inside the photon put to zero. This has a 
dramatic effect. For yi ab < 1, the curves for the direct-photon contribution and the sum 
of the direct and resolved contributions are almost identical. This shows that the resolved 
cross section is essentially made up by the charm content of the photon. Only for yi ab > 1, 
where the cross section decreases with increasing yi a t, can we see the effect of the other 
components of the photon, i.e. the gluon and the light quarks. Therefore charmed-hadron 
photoproduction at large pt seems to be an ideal place to learn specifically about the 
charm content of the photon PJ29[1, which otherwise can only be studied in charmed- 
hadron production in large-Q 2 ej scattering. 

All finite-order perturbative calculations are plagued by scale dependences. Although 
we may expect that this dependence is reduced in NLO, it is not in general negligible, 
in particular at moderate Pt- In Fig. ^, we show the cross section d 2 a j dy^dpT at 
Pt = 7 GeV for our standard PDF choice and scales n = Mf/2 = £m T with £ = 1/2, 1 
and 2. In all these cases, we keep the constraint no = 2m. This scale variation produces 
changes of the cross section of only ±10%, which indicates good perturbative stability. 
The maximum change occurs where the cross section is maximal. The main effect is 
actually due to the variation of the renormalization scale. In fact, if we stick to fi = m?, 
but choose Mf = m^ instead of Mf = 2rriT, we observe only a small difference in the 
cross section; compare the dotted and full curves in Fig. |[ 

The cross section for the production of charmed hadrons depends on the choice of 
FF's. This was already demonstrated in Figs. [l| and 0, where we compared the cross 
sections for different assumptions concerning the FF's. Possibly, the choice e = 0.06 in 
the Peterson FF, which we have adopted from jnj, might not be fully compatible with 
our procedure. In particular, the value of e must depend on the value of the initial scale 
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fiQ. Up to now, we have fixed /j,q = 2m and e = 0.06. In Fig. |7|, we study how variations in 
/i and e affect the cross section d 2 a/dyi a bdpx at px = 7 GeV. If we reduce /i to /i = m 
keeping e = 0.06, the cross section decreases, as may be seen by comparing the dotted 
curve with the full one. This reduction can be compensated by simultaneously adjusting 
e to a smaller value. The variation of the cross section with e may also be inferred from 
Fig. 0, where e = 0.04, 0.06 and 0.08 are considered. Decreasing e leads to a larger 
cross section. This may be understood by observing that case C emerges from case E 
in the limit e — > 0. In order to obtain definite predictions, we must use the information 
on the e parameter in connection with a fixed choice for fio from other D*^ production 
experiments, e.g. from e + e~ — > D*^ + X. So far, the choice e = 0.06 was motivated by 
Chrin's analysis (T^| of charmed-meson production in the PETRA energy range. This 



analysis relied on fragmentation models. Therefore it is unclear whether e = 0.06 is the 
correct choice for the Peterson FF at Ho = 2m. We shall come back to this point when 
we compare our calculation with the experimental data from HI and ZEUS. 

3.4 Comparison with HI and ZEUS data 

In this section, we compare our NLO predictions for the cross section of inclusive D*^ 1 
photoproduction in ep scattering at HERA with three recent sets of data: a) HI tagged, b) 
HI untagged || and c) ZEUS untagged f|. We sum over D* + and D*~ mesons and adjust 
the normalization of the Peterson FF by including the measured branching fraction B(c —>■ 



D* + ) = B{c — > D* ) = 0.260 p7| . We adopt the respective experimental constraints on 
the equivalent photon spectrum of (Q): a) 0.28 < x < 0.65, Q 2 max = 0.01 GeV 2 , b) 
0.1 < x < 0.8, Q 2 max = 4 GeV 2 § and c) 0.147 < x < 0.869, Q 2 max = 4 GeV 2 [|. We first 
consider the px distribution da/dpx integrated over the rapidity interval — 1.5 < y\ a b < 1- 
In Figs. [8]a-c, we confront the HI tagged, HI untagged and ZEUS untagged data with our 
respective NLO predictions. The renormalization and factorization scales are chosen to be 
\i = Mf/2 = £m T with £ = 1/2 (dashed lines), £ = 1 (solid lines) and £ = 2 (dash-dotted 
lines). We observe that the scale dependence is small, approximately ±10%, indicating 
that corrections beyond NLO are likely to be negligible. In all three cases, the agreement 
with the data is remarkably good, even at small px, where we would not have expected it. 
In our approach, the px spectra come out somewhat larger than in the all-massive theory 



1 2(1 which was used for comparison in |||| and therefore agree slightly better with the 



data. In the low-px range, this may be partly attributed to the omission of the power 
terms in m/px in our approach. 

It is well known that the rapidity distribution d 2 a / dyi ab dpx at fixed px usually allows 
for a more stringent test of the theory than da/dpx- Unfortunately, the measured yi a b 
distributions available so far all correspond to d 2 a / 'dyi^dpx integrated over px intervals 
with rather small lower bounds, namely 2.5 GeV < px < 10 GeV in the case of HI and 
3 GeV < px < 12 GeV in the case of ZEUS ||. The corresponding HI tagged, HI untagged 
and ZEUS untagged data are compared with our NLO predictions for £ = 1/2, 1 and 2 in 
Figs. ||a-c, respectively. The theoretical predictions for HI tagged and HI untagged are 
quite different. The curves for the tagged case have their maxima at smaller yi a b, close 
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to yiab = —1.2, and exhibit a much stronger variation with yi^ than the curves for the 
untagged case, which reach their maxima near —0.5. This difference may be understood 
by observing that the soft end of the equivalent photon spectrum is eliminated in the 
tagged case by the cut x > 0.28. The agreement between theory and data is worse than 
in the case of da/dpx- On the other hand, the experimental errors are still rather large. 
We observe that in all three cases there is perfect agreement in the central region of the 
detector, at yi^ = —0.25, whereas in the backward direction the theoretical predictions 
tend to slightly overshoot data, especially in the case of HI. In the case of ZEUS, where 
the lower cut on px is somewhat larger than for HI, our prediction agrees slightly better 
with the data than the result of the all-massive calculation presented in Fig. 4 of ||. 
However, this observation has to be taken with a grain of salt, since the bulk of the 
cross section is accumulated at the low-pr end, where our approach is expected to be less 
reliable. 

In the comparison of the data with our predictions, we have to keep in mind that 
the latter depend on the parameter e and the starting scale /xo of the Peterson FF. We 
have chosen e = 0.06 and /xo = 2m. In Fig. |7|, we have seen that varying e by ±0.02 
changes the cross section by ±15%, which is comparable to the scale dependence. We 
have verified by explicit calculation that the choice e = 0.06 in conjunction with /x = 2m 
yields a satisfactory global fit to the (x) values measured in e + e~ experiments with CM 
energies between 10 and 91.2 GeV. We conclude that, at the present stage, the comparison 
between theory and data is not jeopardized by the uncertainty in e. 

4 Summary and Conclusions 

In this work, we calculated the cross section of inclusive D* ± -meson production via small- 
er 2 ep scattering at HERA energies on the basis of a new massless-charm approach. Specif- 
ically, the MS factorization of the collinear singularities associated with massless charm 
quarks in the final state is adjusted so as to match the corresponding massive-charm cal- 
culation. This is implemented in a way similar to switching from the MS factorization 
scheme e.g. to the deep-inelastic-scattering (DIS) scheme within the massless analysis. 
The fragmentation of the quasi-massive charm quarks into the physical D* + mesons at 
long distances is then described by the NLO-evolved Peterson FF with a suitable choice 
of e at the starting scale /x . Compared to the all- massive calculation, this has the ad- 
vantage that the large logarithms \n(p T /m 2 ) are resummed, in particular through the 
appearance of a much more sizeable resolved-photon contribution. We believe that this 
massless-charm scheme is much more suitable in the large-pr regime. Furthermore, it 
offers us the opportunity to investigate the charm distribution inside the resolved photon 
(see Fig. ||), which does not enter the all- massive calculation. 

We studied various refinements of the massless-charm approach w.r.t. the description 
of the c — > D* + fragmentation. The most realistic fragmentation model, based on the 
Peterson FF evolved through NLO evolution equations from the starting scale /xo = 2m 
up to the characteristic scale of the process, was examined in great detail and used for 
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comparisons with recent data from the HI and ZEUS Collaborations at HERA. We found 
good agreement, in particular for the p? spectra as well as the yi a b spectra in the central 
region. Surprisingly, even at rather low p?, some of our predictions agree with the data 
slightly better than those obtained within the all-massive scheme (compare Figs. |c and 
|9]c with Fig. 4 of 0), which may be an artifact of neglecting the non-singular power terms 
in tti/pt- Since our approach is better justified theoretically for p? ^> m, we hope that 
HERA data on inclusive D* ± photoproduction in the large-pr region will soon be available 
with high statistics. In particular, 7/^ spectra with large minimum-p^ cuts would allow 
for a more stringent test of our fragmentation model in particular and the QCD-improved 
parton model in general, and increase our understanding of the charm distribution inside 
the resolved photon. 

If one determines e from an NLO fit to e + e~ data, the result will depend on whether 
scheme E or scheme F is adopted. According to the factorization theorem, this differ- 
ence is expected to be compensated if one makes NLO predictions for different kinds of 
experiments as long as one works in the respective scheme. By the same token, NLO 
predictions based on the all-massive scheme with Peterson fragmentation [20| suffer from 
an essentially uncontrolled normalization as long as the e parameter is not fitted within 
the same scheme. 

After the completion of this work, new sets of D*^ FF's have been extracted from fits 
to e + e~ data adopting scheme E [30| and a variant of scheme F [31 . 
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A Mellin transform of the Peterson fragmentation 
function 

The Mellin transform of a distribution function, D(x), is defined as 



D(N) = / dxx N - l D(x), 
Jo 

where N is complex. Taking D(x) to be the Peterson FF of ([5|), we find 

2JV 



(6) 



D(N) 



(A+l)(A + 2)(A + 3) 



+2- 



F x (2, 3; N + 4; x^)_ + 2 F X (2, 3; N + 4; x 2 l ) 



x\(x\ - x 2 y 



x\{x x - x 2 ) 2 



2 F 1 (l,3;A + 4;x I 
Xi(Xi - x 2 ) 3 



-1 



2 F 1 {l,3;N + 4;x 2 



x 2 (xi - x 2 ) 3 



(7) 
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where M is the normalization factor in 

Vic) 



2 F l (a, b; c; z) = w^gz h ) J dx ~ x f~ b ~\ l ~ XZ Y 



is the hypergeometric function and 



* V2 = |h±^-i). (9) 

The following expansion is useful for numerical purposes: 

n/ \ r(c)(— z) a m ^ (a) m _)_ n (l — C + a) m +n - n M / \ 

2 F 1 (a, a + m-c-z) = — ' y ' ^ ) ' + z ln(-z 

r(a + m)r(c — a) £f n\[m + n)\ 

+■0(1 + m + n) + -0(1 + n) — ^(a + m + n) — ip(c — a — m — n)] 
a z) -g F (c) ^ r(m-n)(a) B 

+ l j T(a + m) ^ n\T(c-a-n) ' 1 j 

where 

*<*> - < u > 

is the digamma function and 

a „ = ^ ; 12 
r(a) 

is Pochhammer's symbol. 
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FIGURE CAPTIONS 



Figure 1: NLO pr distribution da/dp T , integrated over —1.5 < yi a b < 1, of inclusive D*^ 
photoproduction in ep scattering with untagged positrons as in the ZEUS experiment. 
The branching ratio B(c — > D* + ) is taken to be unity. The fragmentation models A— F 
are described in the text. 

Figure 2: NLO y\ a b distribution d 2 cr/dyi ab dp T , at pt = 7 GeV, of inclusive pho- 
toproduction in ep scattering with untagged positrons as in the ZEUS experiment. The 
branching ratio B(c — > D* + ) is taken to be unity. The fragmentation models A-F are 
described in the text. 

Figure 3: NLO yi a b distribution d 2 <j/dyi a bdpT, at pr = 7 GeV, of inclusive pho- 
toproduction in ep scattering with untagged positrons as in the ZEUS experiment. The 
NLO-evolved Peterson FF is used in connection with the massive-charm factorization 
scheme (case E). The branching ratio B(c — > D* + ) is taken to be unity. The contri- 
butions due to direct (dashed line) and resolved (dash-dotted line) photoproduction are 
shown together with their sum (solid line). 

Figure 4: Influence of the proton PDF's. The total result of Fig. |^ with the CTEQ4M 
set (solid line) is compared with the corresponding calculations with the charm density 
inside the proton switched off (dashed line) and with the MRS(G) set (dash-dotted line). 

Figure 5: Influence of the photon PDF's. The total result of Fig. |3| with the GRV 
set (solid line) is compared with the corresponding calculations with the charm density 
inside the photon switched off (dashed line) and with the GS set (dash-dotted line). For 
comparison, also the direct-photon contribution of Fig. |^ is shown (dotted line). 

Figure 6: Dependence on £, where ft = Mf/2 = (m^. The total result of Fig. |3| with 
£ = 1 (solid line) is compared with the corresponding calculations with £ = 1/2 (dashed 
line) and 2 (dash-dotted line) as well as the one with fi = Mf = itlt (dotted line). 

Figure 7: Dependence on the parameters e and fio of the Peterson FF. The total result 
of Fig. H with e = 0.06 and fio = 2m (solid line) is compared with the corresponding 
calculations with e = 0.04 (dashed line), e = 0.08 (dash-dotted line) and fio = m (dotted 
line). 

Figure 8: The px distributions da/dpx, integrated over —1.5 < yi^ < 1, of inclusive D* ± 
photoproduction in ep scattering with (a) tagged and (b) untagged positrons as measured 
by HI and with (c) untagged positrons as measured by ZEUS are compared with the 
corresponding NLO predictions with £ = 1 (solid lines), £ = 1/2 (dashed lines) and £ = 2 
(dash-dotted lines). 

Figure 9: The yi a b distributions da/dyi a b, integrated over p T , of inclusive D* ± photo- 
production in ep scattering with (a) tagged and (b) untagged positrons as measured by 
HI (2.5 GeV < px < 10 GeV) and with (c) untagged positrons as measured by ZEUS 
(3 GeV < pr < 12 GeV) are compared with the corresponding NLO predictions with 
£ = 1 (solid lines), £ = 1/2 (dashed lines) and £ = 2 (dash-dotted lines). 
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